Photon orbital angular momentum in a plasma vortex 
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We study theoretically the exchange of angular momentum between a photon beam and a plasma vortex, 
and demonstrate the possible excitation of photon angular momentum states in a plasma. This can be relevant 
to laboratory and space plasma diagnostics; radio astronomy self-calibration; and generating photon angular 
momentum beams. A static plasma perturbation with helical structure, and a rotating plasma vortex are studied 
in detail and a comparison between these two cases, and their relevance to the physical nature of photon 0AM, 
is established. 
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It is well known that photons can caiTy not only intrin- 
sic spin angular momentum (SAM), which is associated with 
their polarization state, but also extrinsic orbital angular mo- 
mentum (OAM) UJ. The existence of photon angular mo- 
mentum has always been recognized on theoretical grounds, 
and was first experimentally demonstrated in the 1930's 
While the experimental results led to quite some discussions 
about their proper interpretation fl, they did not excite much 
curiosity at the time as to their utilization. It is also known that 
quantum OAM states are associated with spherical wave func- 
tions |4], and can be excited by pointlike sources. Only re- 
cently photon OAM started receiving considerable attention, 
when it was found that they can be associated not only with 
spherical waves but also with cylindrical waves that can be 
easily produced by laser sources. 

The demonstration that in a laser beam the Laguerre- 
Gaussian modes correspond to well defined OAM modes, and 
that these photon modes cannot only be measured as a pho- 
ton beam property but can also be detected at the single 
photon level [7], explains the relevance of the present research 
on photon OAM states. Utilization of photon OAM states in 
the low frequency (< 1 GHz) radio wave domain, allowing 
digital control of the signals, was recently proposed in Ref.|8l 
as an additional method for characterizing and studying radio 
sources. 

The possibility of studying space plasma vorticity remotely 
by measuring the OAM of radio beams interacting with the 
vortical plasma was pointed out in Ref. |9|. Here we analyze 
this possibility theoretically by studying the exchange of an- 
gular momentum between a plasma medium and a photon 
beam. Such an analysis extends some of the more recent 
studies of optical effects associated with photon OAM into 
the plasma physics domain, and can be relevant for future 
plasma diagnostics, both in laboratory and in space plasma. 
A good understanding of the coupling between plasma vor- 
ticity and radio beam OAM will pave the way for improved 
self-calibration techniques in radio astronomy, and for finding 
new methods of generating electromagnetic beams that carry 



OAM. 

Two plasma situations will be analyzed in detail. First, we 
consider a photon beam propagating in a static plasma pertur- 
bation with a helical vortex structure. We will see that even 
in this simple case, higher OAM states of the photon can be 
excited. Secondly, we consider the interaction of the photon 
beam with a rotating plasma vortex, characterized by a rota- 
tional angular frequency Q.. In contrast with the first case, the 
excitation of higher order OAM states in the second situation 
is accompanied by a frequency shift by multiples of 2Q.. 

In the case of a static plasma vortex, we consider transverse 
electromagnetic waves propagating in a isotropic plasma, 
where ion motions will be neglected. In general terms, we can 
describe these waves by the electric field propagation equation 



(1) 



where the current J — —eny is determined by the electron fluid 
equations 
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-(E + vxB) (2) 



Thermal and relativistic mass effects are ignored. In the pres- 
ence of a static plasma vortex, we can define the mean electron 
density and velocity by 

« = «o + «('", z) cos (Zo(p + q'oz) , v = vo(r,f) + 5v (3) 

where no and vq describe the background plasma conditions 
and 5y is the perturbation associated with the propagating 
electromagnetic wave. We note that the plasma helix vor- 
tex density perturbation, described in cylindrical coordinates 
r = (r, depends on the distance with respect to the vortex 
axis of symmetry and is allowed to vary slowly along z, on a 
scale much longer than the spatial period zq ~ 27r/^o- For a 
typical double vortex, we will have k) = \. Plasma rotation is 
ignored for the moment, but will be considered further below. 
In the case of a static helical perturbation, the plasma cuiTent 
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is J = — eno(r)5v, and the propagation equation can take the 
form 



c2 dfi 



'pO 



1 + e(r, 



E = 



(4) 



where 



CO 



€^fic\ air z) 

e(r,9,z) = ^-^cos(Zo?' + 9oz) (5) 



no 



We further assume wave propagation along the vortex axis Oz, 
and consider solutions of the form 



E(r,f) = A(r)exp -iCOf + ! / k{z)dz 



(6) 



where (O is the wave frequency, and the wave amplitude A(r) 
varies only slowly along z and satisfies 



\d^A/dz^\ « \2kdA/dz\ 



(7) 



We can then reduce the wave equation ^ to the perturbed 
paraxial equation 



cot. 



with 



V] +2ik^-^e{r,(p,z) 
oz 



A = 



(8) 



(9) 



In the absence of the vortex perturbation, equation ([8]) would 
reduce to the usual paraxial optical equation. In this case, a 
general solution can be represented in a basis of orthogonal 
Laguerre-Gaussian modes, according to the expansion 

A(r,(p,z) =^A,,(r,zy''Pexp(-^)e„ (10) 
pi ^ 

where w = w{z) is the beam waist, e^/ are unit polarization 
vectors, and the amphtudes Api are defined by 



Api{r,z) ^Api{z)- 



V2l p\ 
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l^p {x) 



(11) 



where l}'} (x) are the associated Laguerre polynomials, with 
X — /w^. The integers p and I represent the radial and the 
azimuthal (quantum) numbers, respectively. Using equation 
([Tol l in equation (|6]l, we can say that the total electric field is 
represented by a superposition of Laguerre-Gaussian states, in 
the form 



E(r,f) ==52Ep/(r)exp ( -/a)f + / / k{z!)dz: ] (12) 
pi 



with 



(13) 



such that 



Fpiir,(p)°- 



-A ,-,,."^exp 



r 



(14) 



obeying the orthogonality condition 

l-2n 



£rdrj^ d(pF;,Fp,p ^ 5pp,5ip (15) 

In the general case, when a vortex perturbation e{r,(p,z) is 
present, these modes will be coupled to each other, through 
the relation 



p'l' 



(16) 



where the coupling coefficients are defined by 



2k 



K{pI,p'l')^co'poJ^ rdrj^ d(pF;,Fp,,,eir,(p) (17) 

In the simplest case where e is only dependent on the az- 
imuthal angle (p, the coupling coefficients reduce to 



K{pl,p'l') = ajo5, 



'pQOpp' 



2n 



e((p)e'(''-''>1'd(p (18) 



This expression remains valid when the radial scale of the 
plasma vortex is much larger than the photon beam waist w{z)- 
As a special case, we consider a photon beam with no initial 
OAM, which can be described by Ep/ = for / ^ 0. Here the 
mode coupling can be assumed to be sufficiently weak such 
that the zero OAM mode is dominant over the entire inter- 
action region, such that \Epo\ ^ \Ep/ii-^Q\. For a helical static 
plasma perturbation as defined by equation (|4]i, we then obtain 

K{pl,p'l') = ncol^^^dpp, [Szv/o^'"'" + 5/'./o^"*1 (19) 

Replacing this in the coupled mode equation (fTSI l. and inte- 
grating over the axial coordinate z, we obtain, for Ap/(0) = 
A (0)5/0, and assuming the same polarization state for all the 
interacting modes, the following expression for the field mode 
amplitudes 

2c^ Jo k(z') no 

This expression describes the rate of transfer of OAM from the 
static plasma vortex to the electromagnetic field. We notice 
that such a transfer is inhibited after a distance of order In /q^, 
the helical path length. The favorable case is therefore that of 
an interaction distance shorter than this length, or in the limit, 
a plasma structure with no axial periodicity. This picture will 
qualitatively change for the case of a time dependent plasma 
perturbation considered below. 

Equation ( l20l l is only valid when the transfer of OAM is 
small, such that the amplitude of the initial Gaussian mode 
ApQ can be considered constant along the axis. In order to 
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derive a more general solution where the amplitude of the ini- 
tially excited mode is allowed to change, we can go back to 
the coupled mode equations ( fT6] l and ( fTTj l and assume a nearly 
constant coupling coefficient. Writing these equations in a 
simplified but obvious new notation where we drop the radial 
index p, we obtain 



and 



with 



no 



assumed constant along the interaction length. For q{) 
this can be solved in terms of Bessel functions 

Ai,+vi,{z)^r''A{0)Jv{2KQz) 



(21) 



(22) 



0, 



(23) 



where is the initial orbital angular momentum state of the 
electromagnetic beam, and v is an integer This solution 
clearly shows the decay of the initial state over all the other 
states Hi + vlo), on a length scale approximately determined 
by the inverse of the coupling constant Kq. 

Secondly we consider the case of rotating plasma vortex, 
where the plasma density profile is assumed to be homoge- 
neous and constant, e{r,(p) = 0, but the plasma is allowed 
to rotate with respect to the laboratory frame (the frame at- 
tached to the photon beam source), with an angular velocity 
£2, around the axis Oz- This means that plasma particles lo- 
cated at the transverse coordinates (r, (p) have an unperturbed 
velocity given by 

vo = vo(r)ee = vo(r)[-sin(£2f + ^)e.T + cos(£2f + ^)ey] (24) 

In order to determine the perturbed velocity associated with 
the wave field, 8\, we linearize the electron momentum equa- 
tion with respect to this quantity, and write 



— +VO-V 5v + 5v-Vvo = (E + voxB) (25) 

dt I m 



For very slow rotation velocities, we can neglect the magnetic 
field term. This term would lead to a current contribution, 
oscillating at frequencies o ± £1, which is ignored here. The 
remaining equation shows that the plasma rotation introduces 
a frequency shifted term, oscillating at plus or minus twice the 
rotation frequency f2. As a consequence, the resulting plasma 
current has three different frequency components, and can be 
written as 

J(f ) = -enov(f ) = ^ Jv exp[-/(a) + 2vD.)t + 2/ v^] 

v=0,±l 

(26) 

which introduces a coupling between three different fre- 
quency modes o) and ((B±2£2). For a circularly polarized 
electric field E = Zi^ei, we can write 



1 ± e± 

&m V 2(0 dr 



(27) 



l±i 



±i 



e^no 



dvQ 



2V2com (coitn) dr 



■e± 



(28) 



Here we have defined e± = (e^ ± iCy)/ \/2. For arbitrary polar- 
ization, we use E — £+6+ +£'_e_. Notice that the polariza- 
tion states e± correspond to the two photon spin states. The 
above expressions show that a rotating plasma responds dif- 
ferently to the different spin states. More interestingly, these 
spin states become coupled by the frequency shifted terms de- 
termined by J±. Here however, we concentrate our attention 
on the photon OAM states, and consider a linearly polarized 
field, with equal contributions from the two spin states. 

These expressions for the electron current suggest the use 
of the following wave solution 



E(r,f) =I^E„,,,(r,f)exp 

npl 



'i{(0 + nQ.)t + i I k{z')dz' 



with 



(29) 
(30) 



Emp/(r,0 = Kpiiz)Fpi{r, (p) 

where /Vip are determined by equation (fT4l i. This is a straight- 
forward generalization of the solutions (fT2] ) and ( flj] ). Replac- 
ing this form of solution in the modified paraxial equation, 
we obtain new coupled mode equations for the slowly varying 
field amplitudes 



-A , 

dz 



where 



2k„c' 



X I Mnpiyp'l')A„,,,,'e"^''-'-'"^' (31) 



V n'p'l' 

1 

c 



(32) 



The coupling coefficients are now defined as 
Kv{npl, n'p'l') 



2k 



^ -ia„-v J^^ rdr d(p(7y{r)F*p,F,,p,,,e'''"(' (33) 

where the quantities <Tv (r) possibly depend on the radial coor- 
dinate through the rotation frequency £2(r). We assume here 
a constant rotation frequency (over the beam waist distance) 
and also assume linear polarization of the incident wave beam 
in the Ox direction, which leads to 



CO 



Co - 



CO 



CO 



PO 



(34) 



(«±£2) dr 
We are then left with 

Kv{npl,n' p'l') = —2nico„i_vav5ppi5iji+2v (35) 

Replacing this in the coupled mode equations, we get, for a 
well defined but unspecified quantum number p, the following 
result 



-^A„pi — KnvA„j^2ve 

"Z y^^l 



iAvZ 



(36) 
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with Ay = ~ k„), and 



. (On-V 



(37) 



This can be simpUfied further, by noting that the vortex ro- 
tation frequency is much lower than the incident wave fre- 
quency Q.<^ CO, which allows us to replace, for moderate val- 
ues of the integer n, the coupling coefficients ^„v by a constant 
value, K = {ncOpQ/cco^)(dvo/dr), leading to 



■^A{n,l)=K\A{n-2,l+2)e-"'~+A(n+2,l-2)e"^~\ (38) 
oz 

with A ~ 2Q./c. In contrast to the case of static perturbations 
(which only couples states with equal frequency but different 
orbital angular momentum), the case of a rotating but uniform 
plasma couples both frequency and angular momentum states. 
According to equation dSST l. a given mode {(On, I) can decay 
into the modes (a)„+2,' — 2) and (a)„_2,^ + 2). 

For initial photon states with frequency CO = COq and no or- 
bital angular momentum {Iq — 0), we obtain from equation 
I the following cascading decay equations 



dz 



Aa^K 



Aa+ie 



-iAz 



(39) 



with Aa =A{n = 2a, I = —2a). Again, in the limit of a neg- 
ligible dephasing, A ^ 0, we can get simple Bessel function 
solutions, determined by 



A(«o + 2«,/o-2«) = (-l)"A(0)y„(2^z) 



(40) 



where A(0) = A{no,lo) is the amplitude of the initial photon 
state. Notice the similarity with our previous solutions for a 
static vortex. 

Finally, we could consider plasma vortex perturbations with 
some parallel wavelength, of the form exp(/^oz)- In this case, 
we would have to replace A by (A — qojz in the previous mode 
coupled equations. This means that a perfect phase matching 
between the various photon modes can be achieved (at least in 
the limit Q. ^ co„), for A — qo. 

In this work we have considered photon orbital angular mo- 
mentum, for electromagnetic waves propagating in a plasma. 



To our knowledge, this is the first theoretical work dealing 
with photon OAM in the context of plasma vorticity. We have 
studied the propagation of electromagnetic wave beams in an 
isotropic plasma, and we have discussed two distinct physi- 
cal conditions: a helical disturbance in static plasma, and a 
rotating plasma vortex. We have shown that, in both cases, a 
cascading process of OAM transfer between the plasma and a 
photon beam can be achieved. For static plasma perturbations 
with a finite helicity there is no photon frequency shift. Only 
the total angular momentum of the photon beam is modified. 
In contrast, for a rotating plasma perturbation without any he- 
lical perturbation, photon OAM states are excited with photon 
frequency shifts that are multiples of twice the plasma rotation 
frequency. These distinct features can be used as an additional 
diagnostic method of the plasma properties, and can be useful 
in the context of both space and laboratory plasma. 
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